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Abstract. For any finite simple graph G = (V, E), the discrete Dirac operator 
D = d+d* and the Laplace-Beltrami operator L = dd* +d*d = (d+d*) 2 = D 2 
on the exterior algebra bundle f! = ©fi^ are finite v X v matrices, where 
dim(H) = v = ^Zk v k> with V/, = dim(f!i ! ) denoting the cardinality of the set 
of complete subgraphs of G. We prove the McKean-Singer formula 
x(G) = str(e — tL ) which holds for any complex time t, where x(G) = str(l) = 
J^ k (— l) k vj i is the Euler characteristic of G. The super trace of the heat 
kernel interpolates so the Euler-Poincare formula for t = with the Hodge 
theorem in the real limit t — > oo. More generally, for any continuous complex 
valued function / satisfying /(0) = 0, one has the formula x(G) = str(e f ( D '). 
This includes for example the Schrodinger evolutions x(C) = str(cos(tD)) 
on the graph. After stating some immediate general facts about the spectrum 
which includes a perturbation result estimating the Dirac spectral difference of 
graphs, we mention as a combinatorial consequence that the spectrum encodes 
the number of closed paths in the simplex space of a graph. We give a couple 
of worked out examples and see that McKean-Singer allows to find explicit 
pairs of nonisometric graphs which have isospectral Dirac operators. 

1. Introduction 

Some classical results in differential topology or Ricmannian geometry have ana- 
logue statements for finite simple graphs. Examples are Gauss-Bonnet [17] . Poincare- 
Hopf [12], Riemann-Roch |2], Brouwer-Lefschetz jTHj or Lusternik-Schnirelman [15] . 
While the main ideas of the discrete results are the same as in the continuum, there 
is less complexity in the discrete. 

We demonstrate here the process of discretizing manifolds using graphs for the 
McKean-Singer formula [24] 

X(G) = str(e- tL ) , 

where L is the Laplacian on the differential forms and x(C) is the Euler char- 
acteristic of the graph G. The formula becomes in graph theory an elementary 
result about eigenvalues of concrete finite matrices. The content of this article is 
therefore teachable in a linear algebra course. Chunks of functional analysis like 
elliptic regularity which are needed in the continuum to define the eigenvalues are 
absent, the existence of solutions to discrete partial differential equations is trivial 
and no worries about smoothness or convergence are necessary. All we do here is 
to look at eigenvalues of a matrix D defined by a finite graph G. The discrete 
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approach not only recovers results about the heat flow, we get results about uni- 
tary evolutions which for manifolds would require need analytic continuation: also 
the super trace of U(t) = e ltL where L is the Laplace-Beltrami operator is the 
Euler characteristic. While U(t)f solves (id/dt — L)f = 0, the Dirac wave equation 
(d 2 /dt 2 + L) = (d/dt - W){d/dt + iD)f = is solved by 

(1) f(t) = U(t)f(0) = ax(Dt)f(0) + sin(Dt)D- x f(0) , 

where /'(0) is on the orthocomplcmcnt of the zero eigenspace. Writing ip = 
f — iD^ 1 /' we have ip(t) = e lDt tp(0). The initial position and velocity of the 
real flow are naturally encoded in the complex wave function. That works also for 
a discrete time Schrodinger evolution like T(f(t),f(t-1)) = (Df(t) - f(t~ 1), /(*)) 
[15] if D is scaled. 

Despite the fact that the Dirac operator D for a general simple graph G — (V, E) 
is a natural object which encodes much of the geometry of the graph, it seems not 
have been studied in such an elementary setup. Operators in HU [25J 0] have no 
relation to the Dirac matrices D studied here. The operator D 2 on has appeared 
in a more general setting: [23] builds on work of Dodziuk and Patodi and studies 
the combinatorial Laplacian (d + d*) 2 acting on Cech cochains. This is a bit more 
general: if the cover of the graph consists of unit balls, then the Cech cohomology 
by definition agrees with the graph cohomology and the Cech Dirac operator agrees 
with D. Notice also that for the Laplacian Lq on scalar functions, the factorization 
Lq = rfo^o i s very well known since do is the incidence matrix of the graph. The 
matrices dk have been used by Poincare already in 1900. 

One impediment to carry over geometric results from the continuum to the discrete 
appears the lack of a Hodge dual for a general simple graph and the absence of 
symmetry like Poincare duality. Does some geometric conditions have to be im- 
posed on the graph like that the unit spheres are sphere-like graphs to bring over 
results from compact Riemannian manifolds to the discrete? The answer is no. The 
McKean-Singer supersymmetry for eigenvalues holds in general for any finite simple 
graph. While it is straightforward to implement the Dirac operator D of a graph in 
a computer, setting up D can be tedious when done by hand because for a complete 
graph of n nodes, D is already a (2™ — 1) X (2" — 1) matrix. For small networks 
with a dozen of nodes already, the Dirac operator acts on a vector space on vector 
spaces having dimensions going in to thousands. Therefore, even before we had the 
computer routines were in place which produce the Dirac operator from a general 
graph, computer algebra software was necessary to experiment with the relatively 
large matrices first cobbled together by hand. Having routines which encode the 
Dirac operator for a network is useful in many ways, like for computing the coho- 
mology groups. Such code helped us also to write [13] . Homotopy deformations of 
graphs or nerve graphs defined by a Cech cover allow to reduce the complexity of 
the computation of the cohomology groups. 

Some elementary ideas of noncommutative geometry can be explained well in this 
framework. Let H = L 2 (Q) denote the Hilbert space defined by the simplex set 
Q of the graph and let B{H) denote the Banach algebra of operators on H . It 
is a Hilbert space itself with the inner product (A,B) — tr(A*B). The operator 
D € B(H) defines together with a subalgebra A of B(H), a Connes distance on 
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Q. Classical geometry is when A = C(Q) is the algebra of diagonal matrices. In this 
case, the Connes metric 8(x, y) = sup ;ey!l ^ D ;]|<i(i, s. x — e y ) between two vertices is 
the geodesic shortest distance but S extends to a metric on all simplices Q. The ad- 
vantage of the set-up is that other algebras A define other metrics on Q and more 
generally on states, elements in the unit sphere X oi A* . In the classical case, 
when A = H = C(Q) — R v , all states are pure states, unit vectors v in H and 
correspond to states I — > tr(l ■ e v ) — (lv,v) in A*. For a noncommutative algebra 
A, elements in X are in general mixed states, to express it in quantum mechanics 
jargon. We do not make use of this here but note it as additional motivation to 
study the operator D in graph theory. 

The spectral theory of Jacobi matrices and Schrodinger operators show that the dis- 
crete and continuous theories are often very close but that it is not always straight- 
forward to translate from the continuum to the discrete. The key to get the same 
results as in the continuum is to work with the right definitions. For translating 
between compact Riemannian manifolds and finite simple graphs, the Dirac oper- 
ator can serve as a link, as we see below. While the result in this paper is purely 
mathematical and just restates a not a priory obvious symmetry known already 
in the continuum, the story can be interesting for teaching linear algebra or illus- 
trating some ideas in physics. On the didactic side, the topic allows to underline 
one of the many ideas in [21], using linear algebra tools only. On the physics side, 
it illustrates classical and relativistic quantum dynamics in a simple setup. Only 
a couple of lines in a modern computer algebra system were needed to realize the 
Dirac operator of a general graph as a concrete finite matrix and to find the quan- 
tum evolution in equation ([!]) as well as concrete examples of arbitrary members of 
cohomology classes by Hodge theory. The later is also here just a remark in linear 
algebra as indicated in an appendix. 

2. Dirac operator 

For a finite simple graph G — (V,E), the exterior bundle fl = (Bk^k is a finite 
dimensional vector space of dimension n = X)/c=o Wfc > wnere v k — dim(fik) is the 
cardinality of Qk the set of Kk+\ simplices contained in G. As in the continuum, 
we have a super commutative multiplication A on fi but this algebra structure is 
not used here. The vector space Qk consists of all functions on Qk which are an- 
tisymmetric in all k + 1 arguments. The bundle splits into an even Qf, ~ (Bk^2k 
and an odd part Qf — (Bk^2k+i which are traditionally called the bosonic and 
fcrmionic subspaces of Q. The exterior derivative d : ft — > f2, df(xo, . . . , x n ) — 
J2k=o i~^) k f( x o> ■ ■ ■ ,ik, ■ ■ ■ , £ n -i) satisfies d 2 = 0. While a concrete implementa- 
tion of d requires to give an orientation of each element in Q, all quantities we are 
interested in do not depend on this orientation however. It is the usual choice of a 
basis ambiguity as it is custom in linear algebra. 

Definition 1. Given a graph G = (V, E), define the Dirac operator D = d + d* 
and Laplace-Beltrami operator L = D 2 = dd* + d*d. The operators L p = 
d*d p + dp-id* _ 1 leaving £l p invariant are called the Laplace-Beltrami operators L p 
on p- forms. 

Examples. We write to indicate multiplicity m. 

1) For the complete graph K n the spectrum of D is {~\/n , 0, — ^/n }. 



4 



OLIVER KNILL 




v -30 v - 159 v - 222 v - 95 v - 16 v - 1 



dim= 2 84 v= 523 



bfl= 1 bp 3 b2=15 bj=0 b 4 =0 bj= 

























: 1 '. ■. \ ■ \ 








— : • ■ ■ 












% 








; 























Figure 1 . The lower part shows the Dirac operator D = d + d* 
and the Hodge Laplace-Beltrami operator L = dd* + d*d of a ran- 
dom graph with 30 vertices and 159 edges. The Dirac operator D 
of this graph is a 523 x 523 matrix. The Laplacian decomposes into 
blocks L p , which are the Laplacians on p forms. The eigenvalues of 
L/c shown in the upper right part of the figure are grouped for each 
p. The multiplicity of the eigenvalue is by the Hodge theorem the 
p'th Betti number b p . In this case bo = 1, b\ = 3, 62 = 13. Nonzero 
eigenvalues come in pairs: each fermionic eigenvalue matches a 
bosonic eigenvalue. This super symmetry discovered by McKean 
and Singer shows that str(e itL ) = £ A (_i)dcg(A) e itA = = 
X) p (— l) p b p = X) P (~ l) Py p- While all this is here just linear alge- 
bra, the McKean formula mixes spectral theory, cohomology and 
combinatorics and it applies for any finite simple graph. 

2) For the circle graph C n the spectrum is {±y/2 — 2 cos(2n k/n)}, k = 0, n — 1 , 
where the notation understands that has multiplicity 2. The product of the 
nonzero eigenvalues, a measure for the complexity of the graph is n 2 . 

3) For the star graph S n the eigenvalues of D are {— y/n, — l( n_1 ), 0, l^ 1-1 ), xfn }. 
The product of the nonzero eigenvalues is n. 

4) For the linear graph L n with n vertices and n — 1 edges, the eigenvalues of D are 
the union of {0}, and ±aK, where K is the (n — 1) x (n— 1) matrix which has 2 
in the diagonal and 1 in the side diagonal. The product of the nonzero eigenvalues 
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is — n. 
Remarks. 

1) The scalar part Lq is one of the two most commonly used Laplacians on graphs 
[6]. The operator L generalizes it to all p-forms, as in the continuum. 

2) When dk is implemented as a matrix, it is called a signed incidence matrix of 
the graph. Since it depends on the choice of orientation for each simplex, also the 
Dirac operator D depends on this choice of basis. Both L and the eigenvalues of D 
do not depend on it. 

3) The Laplacian for graphs has been introduced by Kirkhoff in 1847 while proving 
the matrix-tree theorem. The computation with incidence matrices is as old as al- 
gebraic topology. Poincare used the incidence matrix d already in 1900 to compute 
Betti numbers PI38]. 

4) While D exchanges bosonic and fermionic spaces D : — > rif,ttf — > il^, the 
Laplacian L is the direct sum Zfc : £7^ — > Qk of Laplacians on k forms. 

5) The kernel of Lf~ is called the vector space of harmonic k forms. Such har- 
monic forms represent cohomology classes. By Hodge theory (see appendix) the 
dimension of the kernel is equal to the fc'th Betti number, the dimension of the fc'th 
cohomology group H k {G) = ker(dk)/im(dk-i). 

6) Especially, Lq is the graph Laplacian B — A, where B is the diagonal matrix 
containing the vertex degrees and A is the adjacency matrix of the graph. The 
matrix Lq is one of the natural Laplacians on graphs |6] . 

7) While d p _idp_ x always has p + 1 in the diagonal, the diagonal entries d*d p (x, x) 
counts the number of K p+ \ simplices attached to the p simplex x and d p d*_ 1 (x, x) 
is always equal to p + 1. The diagonal entries Lk(x,x) therefore determine the 
number of k + 1 dimensional simplices attached to a A; dimensional simplex x. We 
will see that the Dirac and Laplacian matrices are useful for combinatorics when 
counting closed curves in Q in the same way as adjacency matrices are useful in 
counting paths in G. 

8) The operator D can be generalized as in the continuum to to bundles. The 
simplest example is to take a 1-form A G SI 1 and to define cZa/ = df + A A /. 
This is possible since 17 is an algebra. We have then the generalized Dirac operator 
Da = d>A + d* A and La — D A . The curvature operator Fg = 4° is no more 
zero in general. McKean-Singer could generalizes to this more general setup as Da 
provides the super-symmetry. 

9) A 1-form A defines a field F = dA which satisfies the Maxwell equations 

dF = 0, d* F = j . Physics asks to find the field F, given j. While this corresponds 
to the 4-current in classical electro magnetism which includes charge and electric 
currents, it should be noted that no geometric structure is assumed for the graph. 
The Maxwell equations hold for any finite simple graph. It defines the evolution of 
light on the graph. The equation d*dA = j can in a Coulomb gauge where d* A = 
be written as L±A = j where L\ is the Laplacian on 1 forms. Given a current j, 
we can get A and so the electromagnetic field F. by solving a system of linear 
equations. This is possible if j is perpendicular to the kernel of L\ which by Hodge 
theory works if G is simply connected. Linear algebra determines then determines 
the field F, a function on triangles of a graph. As on a simply connected compact 
manifold, a simply connected graph does not feature light in vacuum since LA = 
has only the solution A = 0. 
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10) Dirac introduce the Dirac operator D in the continuum to make quantum me- 
chanics Lorentz invariant. The quantum evolution becomes wave mechanics. Al- 
ready in the continuum, there is mathematically almost no difference between the 
wave evolution with given initial position and velocity of the wave the Schrodingcr 
evolution with a complex wave because both real and imaginary parts of e lDt (u — 
iD~ 1 v) = cos(Dt)u + sin(Dt)D~ 1 v solve the wave equation. The only difference 
between Schrodingcr and wave evolution is that the initial velocity v of the wave 
must be in the orthocomplement of the zero eigenvalue. (This has to be the case 
also in the continuum, if we hit a string, the initial velocity has to have average 
zero momentum in order not to displace the center of mass of the string.) A similar 
equivalence between Schrodinger and wave equation holds for a time-discretized 
evolution (u, v) — > (v — Du, u). 

11) An other way to implement a bundle formalism is to take a unitary group 
U(N) and replace entries 1 in d with unitary elements U and —1 with — U. Now 
D = d + d* and L = D 2 are sclfadjoint operators on L 2 (Q,C N ) for some N and 
can be implemented as concrete v ■ N matrices. The gauge field U defines now 
curvatures, which is a J7(iV)-valued function on all triangles of Q. Again, the more 
general operator D is symmetric and supersymmetry holds for the eigenvalues. Ac- 
tually, as tr(D n ) can be expressed as a sum over closed paths, the eigenvalues of Da 
are just N copies of the eigenvalues of D if the multiplicative curvatures are 1 (zero 
curvature) and the graph is contractible. Since tr(/(D)) = for odd functions / 
and tr(D 2 ) is independent of U, the simplest functional on U which involves the 
curvature is the Wilson action tr(Z) 4 ) which is zero if the curvature is zero. It is 
natural to ask to minimize this over the compact space of all fields A. More general 
functionals are interesting like det*(L), the product of the nonzero eigenvalues of 
L; this is an interesting integral quantity in the flat U = 1 case already: it is a 
measure for complexity of the graph and combinatorially interesting since det*(Lo) 
is equal to the number of spanning trees in the graph G. 



3. McKean-Singer 

Let G = (V, E) be a finite simple graph. If K n denotes the complete graph with n 
vertices, we write Q n for the set of K n+1 subgraphs of G. The set of all simplices 
Q = U n G n is the super graph of G on which the Dirac operator lives. The graph 
G defines Q without addition of more structure but the additional structure is 
useful, similarly as tangent bundles are useful for manifolds. It is useful to think of 
elements in Q as elementary units. If v n is the cardinality of Q n then the finite sum 

X (G) = jr(-i) n v n 

n=0 

is called the Euler characteristic of G. 

The set il n of antisymmetric functions / on Q n is a linear space of fc-forms. The 
exterior derivative d : VL n — > £l n +i defined by 



df{x , xx,..., x n+1 ) = ^(-l) fe /(a;o, ...,x k ,..., x n+1 ) 
k 
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is a linear map which satisfies d o d = 0. It defines the cohomology groups 
H n (G) — ker(d n )/im(o? n _ 1 ) of dimension b n (G) for which the Euler-Poincare for- 
mula x(G) = J2 n (~ l) n bn{G) holds. The direct sum SI = © n f2„ is the discrete 
analogue of the exterior bundle. 

McKean and Singer have noticed in the Riemannian setup that the heat flow e~ tL 
solving (d/dt + L)f = has constant super trace str(L(t)). The reason is that 
the nonzero bosonic eigenvalues can be paired bijectively with nonzero fermionic 
eigenvalues. It can be rephrased that D is an isomorphism from Sl£ — > Sl^ , where 
Sl^ is the orthogonal complement of the zero eigenspace in Sib and Sl~^ the orthog- 
onal complement of the zero eigenspace in Slf. In the limit t —> oo we get the 
dimensions of the harmonic forms; in the limit t — > we obtain the simplicial Eulcr 
characteristic. 

Theorem 1 (McKean-Singer). For any complex valued continuous junction f on 
the real line satisfying /(0) = we have 

X (G) = str(cxp(f(D))) . 

Especially, str(exp(— tL)) = x(G) for any complex t. 

Proof. Since D = d + d* is a symmetric v x v matrix, the kernel of D and the kernel 
of L = D 2 are the same. The functional calculus defines f{D) for any complex 
valued continuous function. Since D is normal, DD* = D*D = —L, it can be 
diagonalized D = U*ED with a diagonal matrix E and defines f(D) = U*f(E)U 
for any continuous function. Because / can by Weierstrass be approximated by 
polynomials and the diagonal entries of D 2k+1 are empty, we can assume that 
f(D) = g(D 2 ) = g(L) for an even function g. Let il + be the subspace of SI spanned 
by nonzero eigenvalues of L. Let f2+ be the subspace generated by eigenvectors to 
nonzero eigenvalues of L on Sl p and f2j be the subspace generated by eigenvectors 
to nonzero eigenvalues of L on Slf. Since D commutes with L, each eigenvector 
/ of g(L) on S}f has an eigenvector Df of g{L) on Sl p . Since D : fi+ — > fi+ is 
invertible, there is a bijection between fermionic and bosonic eigenvalues. Each 
nonzero eigenvalue appears the same number of times on the fermionic and bosonic 
part. □ 

Remarks. 

1 By definition, str(l) = X)k=o(~l) kv k agrees with the Euler characteristic. It can 
be seen as an analytic index ind(D) of the restricted Dirac operator Sib — > tlf 
because ker(D) is the space of harmonic states in Sib and coker(D) = ker(D*) is the 
dimension of the fermionic harmonic space. 

2) The original McKean-Singer proof works in the graph theoretical setup too as 
shown in the Appendix. 

3) In |15j we have for numerical purposes discretized the Schroedinger flow. This 
shows that we can replace the flow e tD by a map T(f, g) = (g — Df, f) with 
a suitably rescaled D which is dynamically similar to the unitary evolution and 
has the property that the system has finite propagation speed. The operator 
T 2 (f,g) = (/- D(g -Df),g- Df) = (/ + Lf,g) - (Dg,Df) has the super trace 
str(T 2 ) = str((l, 1)) = 2\(G) so that also this discrete time evolution satisfies the 
McKean-Singer formula. 

4) The heat flow proof interpolating between the identity and the projection onto 
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harmonic forms makes the connection between Euler-Poincare and Hodge more 
natural. While for t — and t — oo we have a Lefshetz fixed point theorem, for 
< t < oo it can be seen as an application of the Atiyah-Bott generalization of 
that fixed point theorem. In the discrete, Atiyah-Bott is very similar to Brouwer- 
Lefshetz [18]. For t = the Lefshetz fixed point theorem expresses the Lefshetz 
number of the identity as the sum of the fixed points by Poincare-Hopf. In the 
limit t — > oo, the Lefshetz fixed point theorem sees the Lefshetz number as the 
signed sum over all fixed points which are harmonic forms. For positive finite t we 
can rephrase McKean-Singer's result that the Lefshetz number of the Dirac bundle 
automorphism eT tL is time independent. 

4. The spectrum of D 
We start with a few elementary facts about the operator D: 

Proposition 2. Let X denote the eigenvalues of D and let deg denote the maximal 
degree of G. If X is an eigenvalue, then —X is an eigenvalue too so that E[A] = 
E<Ai = 0. 



Proof. If we split the Hilbert space as f2 = fi p © f2/ then D = 



A* 
A 



where 



A = d + d* is the annihilation operator and A* which maps bosonic to fermionic 

states and A is the creation operator. Define P = J ^ Then L = Z) 2 , P 2 = 

1,DP + PD = is called supersymmetry in space dimensions (see [3^1 [7]). If 
Df = Xf, then PDPf = -Df = -A/. Apply P again to this identity to get 
D(Pf) — —X(Pf). This shows that Pf is an other eigenvector. □ 

Proposition 3. Every eigenvalue X is contained in [— \/2deg, -\/2degJ so that ev- 
ery eigenvalue of L = D 2 is contained in the interval [0, 2deg] and Var[A] = 
£lLiA 2 /n<2deg. 

Proof. Because all entries of D are either —1 and 1 and because there are d such 
entries, each row or column has maximal length \fd. That the standard deviation is 
< i/2deg follows since the "random variables" Xj take values in [— v / 2dcg, \/2deg]. 
Also a theorem of Schur confirms that J^i — J2i j \^ij\ 2 — n ' 2deg. □ 

Remarks. 

1) It follows that the spectrum of D is determined by the spectrum of L. If 
Hj € [0, deg] are the eigenvalues of L, then ±^/m,j are the eigenvalues of D. It fol- 
lows already from the reflection symmetry of the eigenvalues of D that the positive 
eigenvalues of L appear in pairs. The McKean-Singer statement is stronger than 
that. It tells that if one member of the pair is in the bosonic sector, the other is in 
the fermionic one. 

2) The Schur argument given in the proof is not quite irrelevant when looking at 
spectral statistics. Since the average degree in the Erdoes-Renyi probability space 
E(n, 1/2) is of the order log(n) the standard deviation of the eigenvalues is of the 
order 21og(n) even so we can have eigenvalues arbitrarily close to y/2(n — 1). 

3) Numerical computations of the spectrum of D for large random matrices is dif- 
ficult because the Dirac matrices are much larger than the adjacency matrices of 
the graph. It would be interesting however to know more about the distribution of 
the eigenvalues of D for large matrices. 
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4 Sometimes, an other Laplacian Kq defined as follows for graphs. Let V$(x) 
denote the degree of a vertex x. Define operator K xx = 1 if Vq(x) > and 
K xy — — (Vo(x)Vo(j/)) 1 ^ 2 if ( x ,y) is an edge. This operator satisfies K = CLqC, 
where C is the diagonal operator for which the only nonzero entries are Vq(x)^ 1 ^ 2 
if V (x) > 0. While the spectrum of L is in [0, 2deg] the spectrum of K is in [0, 2]. 
Since Lq has integer entries, it is better suited for combinatorics. In the following 
examples of spectra we use the notation indicating that the eigenvalue A ap- 
pears with multiplicity n. 

For the complete graph K n+ \ the spectrum of Kq is {0, (n + l)/n)( n '} while 
the spectrum of L is {0, }. For a cycle graph C n the spectrum of K is 
{1 — cos(27rfc/n) } while the spectrum of L is 2 — 2 cos(2nk/n). This square graph 
C4 shows that the estimate a(L) C [0, 2deg] is optimal. For a star graph S n which is 
an example of a not a vertex regular graph, the eigenvalues of K are {0, l(™~ 2 ', 2 } 
while the spectrum of L is {0, \( n ~ 2 \n }. 



Proposition 4. The number of zero eigenvalues of D is equal to the sum of all the 

Betti numbers b k . 

Proof. This follows from the Hodge theorem (see Appendix) : the dimension of the 
kernel of L k is equal to b k . □ 

Definition 2. The Dirac complexity of a finite graph is defined as the product 
of the nonzero eigenvalues of its Dirac operator D. 

The Euler-Poincare identity assures that J2iLo(~^-Y( Vi ~ &*) = 0- ^ wc ig nore the 
signs, we get a number of interest: 

Lemma 5. The number of nonzero eigenvalue pairs in D is the sign-less Euler- 
Poincare number J2'iLo( Vi ~~ bi)/2. 

Proof. The sum ^ i Vi — v is the total number of eigenvalues and by the previous 
proposition, the sum J2i h is the total number of zero eigenvalues. Since the eigen- 
values of D come in pairs ±A, the number of pairs is the sign-less Euler-Poincare 
number. □ 

Corollary 6. The sign-less Euler-Poincare number is even if and only if the Dirac 
complexity is positive. 

Proof. Arrange the product of nonzero eigenvalues of D as a product of pairs — X 2 . 

□ 

Corollary 7. If G is a triangularization of a sphere and which has an even number 
of edges, then the Dirac complexity is positive. 

Proof. We have x(G) = v — v 1 + v 2 = 2 and since 6 = b 2 = 1 and b\ = we have 
bi = 2. We can express now the sign-less Euler-Poincare number in terms of the 
number of edges: 

[(«o + v t + v 2 ) - (b + 61 + b 2 )]/2 = [(2 + 2 Vl ) - 2]/2 = 2^/2 = v x . 

□ 

The following explains why the dodecahedron or cube have negative complexity: 
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Corollary 8. If G is a connected graph without triangles which has an even number 
of vertices, then the Dirac complexity is negative. 

Proof. x{G) — Vo — v\ = bo — b\ = 1 — b\. The sign-less Euler-Poincare number is 
[(«o+«i)-(6o + 6i)]/2 = [l-6i+2«i-(l + 6i)]/2 = t;i-6 1 = «o-6o=«o-l. □ 

For cyclic graphs C n the complexity is n 2 if n is odd and — n 2 if n is even. For star 
graphs S n , the complexity is n is odd and — n if n is even. 

Corollary 9. For a tree, the Dirac complexity is positive if and only if the number 
of edges is even. 

Proof. The Dirac complexity is still v\ — b\ as in the previous proof but now b\ = 
so that it is v\. □ 

We have computed the complexity for all Platonic, Archimcdian and Catalan solids 
in the example section. All these 31 graphs have an even number v\ of edges and 
an even number v\ of vertices. 



5. Perturbation of graphs 

Next we estimate the distance between the spectra of different graphs: We have 
the following variant of Lidskii's theorem which I learned from 22 : 

Lemma 10 (Lidskii). For any two self adjoint complex n x n matrices A and B 
with eigenvalues ot\ < ct2 < • • ■ < ot n and fi\ < f$i < • • • < /?„, one has 

n n 

EK-AI^ E A - B \n ■ 
j=i j,i=i 

Proof. Denote with 7, £ R the eigenvalues of the selfadjoint matrix C := A — B 
and let U be the unitary matrix diagonalizing C so that Diag(7i, . . . , j n ) = U CU*. 
We calculate 

i i i,k,l 

< E ■ lEc- 1 ) 7 " 4 ^^! < Ei^'i • 

k,l i k,l 

The claim follows now from Lidskii's inequality J^j \ a j~ Pj\ ^ Ylj ( see ES]) ^ 

This allows to compare the spectra of Laplacians Lq of graphs which are close. Lets 
define the following metric on the Erdoes-Renyi space G[n) of graphs of order n on 
the same vertex set. Denote by d (G,H) the number of edges at which G and H 
differ divided by n. Define also a metric between their adjacency spectra A, ft as 

1 " 

d (A,/2) = -E' A J ■ 
" i=i 

Corollary 11. If the maximal degree in either G,H is deg, then the adjacency 
spectra distance can be estimated by 

d (A,M)=2degdo((7o(G),ao(H)) . 
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Proof. We only need to verify the result for do(G, H) = 1, since the left hand side 
is the I 1 distance between the vectors A, ^ and the triangle inequality gives then 
the result for all d (G, H). For d (G, H) = 1, the matrix C = B - A differs by 1 
only in 2deg entries. The Lidskii lemma implies the result. □ 

Remarks. 

1) We take a normalized distance between graphs because this is better suited for 
taking graph limits. 

2) As far as we know, Lidskii's theorem has not been used yet in the spectral theory 
of graphs. We feel that it could have more potential, especially when looking at 
random graph settings Hi . 

3) For the Laplace operator L Q this estimate would be more complicated since the 
diagonal entries can differ by more than 1. It is more natural therefore to look at 
the Dirac matrices for which the entries are only 0, 1 or —1. 

To carry this to Dirac matrices, there are two things to consider. First, the matrices 
depend on a choice of orientation of the simplices which requires to chose the same 
orientation if both graphs contain the same simplex. Second, the Dirac matrices 
have different size because D is a v x v matrix if v is the total number of simplices 
in G. Also this is no impediment: take the union of all simplices which occur 
for G and H and let v denote its cardinality. We can now write down possibly 
augmented v x v matrices D(G), D(H) which have the same nonzero eigenvalues 
than the original Dirac operator. Indeed, the later matrices are obtained from 
the augmented matrices by deleting the zero rows and columns corresponding to 
simplices which are not present. 

Definition 3. Define the spectral distance between G and H as 

1 

where \j,(J-j are the eigenvalues of the augmented Dirac matrices G,H, which are 
now both v x v matrices. 

Definition 4. Define the simplex distance d(G, H) of two graphs G, H with 
vertex set V as (1/v) times the number of simplices of G,H which are different 
in the complete graph on V . Here v is the total number of simplices in the union 
when both graphs are considered subgraphs of the complete graph on V . 

Definition 5. The maximal simplex degree of a simplex x of dimension k in a 
graph G is the sum of the number of simplices of dimension k + 1 which contain x 
and the sum of the number of simplices of dimension k — 1 which are included in x. 

In other words, the maximal simplex degree is the number of nonzero entries d XiV 
in a column x of the incidence matrix d. We get the following perturbation result: 

Corollary 12. The Dirac spectra of two graphs G,H with vertex set V satisfies 

d(a(G),a(H)) < 2deg ■ d(G,R) , 
if deg is the maximal simplex degree. 

Proof. Since we have chosen the same orientation for simplices which are present 
in both graphs, this assures that the matrix D{G) — D(H) has entries of absolute 
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value < 1 and are nonzero only at D x ^ yi where an incidence happens exactly at one 
of the graphs G, H. The value 2degd(G, H) is an upper bound for - J2k i \D{G)ki — 
D(H)ki\, which by Lidskii is an upper bound for the spectral distance ^ J2^=i IA7 — 
Mil- ' "'" □ 

Remarks. 

1) While the estimate is rough in general, it can become useful when looking at 
graph limits or estimating spectral distances between Laplace eigenvalues of differ- 
ent regions. If we look at graphs of fixed dimensions like classes of triangularizations 
of a manifold M, then deg is related to the dimension of M only and both sides of 
the inequality have a chance to behave nicely in the continuum limit. 

2) For large graphs which agree in many places, most eigenvalues of the Laplacian 
must be close. The result is convenient to estimate the spectral distance between 
the complete graph K n and G. In that case, 2 • deg = 2(2 n_1 — 1) < 2 n and we have 

d(a(G),a(K n ))<d(G,K n ) . 

Since the Dirac spectrum of K n is contained is the set {— y/n, 0, y/n }. It follows that 
Erdoes-Renyi graphs in G(n,p) for probabilities p close to 1 have a Dirac spectrum 
concentrated near ±y/n. If G has m simplices less than K n , then d(G, K n ) < m/2 n . 

Examples. 

1) Let G be the triangle and H be the line graph with three vertices. Then v = 7. 
We have deg = 3. Since the graphs differ on the triangle and one edge only, we 
have d{G, H) = 2/7 and 2dcg • d(G, H) = 12/7 = 1.71 .... The augmented Dirac 
matrices are 
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The spectrum of the Dirac matrix of G is A — {— \/3, — v/3, — \/3, 0, y/3, V3, V3 } 
The spectrum of the augmented Dirac matrix of H is \i = {—■ v/3, — 1, 0, 0, 0, 1, V3 }. 
The spectral distance is d(a(G),a(H)) = (2 + 2\/S)/7 = 0.781 .... 

2) If G is K n and H is K n _ x then v = 2 n - 1. We have d(G,H) = [(2" - 1) - 
(2"- 1 - l)]/(2™ - 1) = 2""7(2" - 1) - 1/2 and 2deg = (2(2 11 - 1 - 1)) ~ 2" and 
so 2degd(G,H) <~ 2^n — 1). The Dirac eigenvalues differ by \^/n — \/n—l)\ at 
(2™- 1 - 2) places and by y/n at (2 n - 1) - (2™- 1 - 1) = 2™- 1 places. The spectral 

distance is \y/n — yn — 1)| 2 "-i V"- 2"-i wn i cn is about yn/2. Taking away 
one vertex together with all the edges is quite a drastic perturbation step. It gets 
rid of a lot of simplices and changes the dimension of the graph. 

3) Let G be the wheel graph W4 with 4 spikes. It is the simplest model for a planar 
region with boundary. Let H be the graph where we make a pyramid extension over 
one of the boundary edges. This is a homotopy deformation. The graph G has 17 
simplices and the graph H has v = 21 simplices. The graph H has one vertex, two 
edges and one triangle more than G so that d(G, H) — 4/21 = 0.190 .... The maxi- 
mal degree is deg = 4 so that the right hand side of the estimate is 32/21 = 1.52 
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The eigenvalues of the Dirac operator D(H) = 



(where, wc wrote a = — 1 for typographic reasons) of H are {—2.370, —2.302, —2.266, — 
, -1.913, -1.839, -1.732, -1.5, -1.302, -0.9296, 0, 0.9296, 1.302, 1.529, 1.732, 1.839, 1.913, 
2,2.266,2.302,2.370}. The eigenvalues of the augmented Dirac operator of G are 
{-2.236,-2.236,-2.236 ,-1.732, -1.732,-1.732, -1.732,-1,0,0, 0,0,0,1 ,1.732,1.732, 
1.732, 1.732, 2.236, 2.236, 2.236}. The spectral differences is 0.337998 .... We were 
generous with the degree estimate. The new added point adds simplices with max- 
imal degree 3 so that the proof of the perturbation result could be improved to get 
an upper bound 8/7 = 1.142 .... It is still by a factor 3 larger than the spectral 
difference, but we get the idea: if H is a triangularization of a large domain and we 
just move the boundary a bit by adding a triangle then the spectrum almost does 
not budge. This will allow us to study spectral differences lim„^oo j- X^j=i I A? ~ Mj I 
of planar regions in terms of the area of the symmetric difference. 

6. Combinatorics 

Adjacency matrices have always served as an algebraic bridge to study graphs. The 
entry A™ has the interpretation as the number of paths in G starting at a vertex i 
and ending at a vertex j. From the adjacency matrix A, the Laplacian Lq = B — A 
is defined. Similarly as for Lq, we can read off geometric quantities also from the 
diagonal entries of the full Laplacians L = D 2 on p forms. 

Definition 6. The number of (p + l)-dimensional simplices which contain the p- 
dimensional simplex x is called the p-degree of x. It is denoted deg p (x). 

Examples. 

1) For p = 0, the degree deg (x) is the usual degree vertex degree deg(x) of the 
vertex x. 

2) For p = 1 the degree deg 1 (x) is the number of triangles which are attached to 
an edge x. 

Proposition 13 (Degree formulas). Forp > we have deg p (x) = L p (x,x) — (p+l). 
For p — we have deg (a;) = L (x,x). 

Proof. Because d p has p + 1 nonzero entries 1 or —1 in each row x, we have 
d*d p (x, x) = p + 1. The reason for the special situation p — is that L is the only 
Laplacian which does not contain a second part dod^ because d^ is zero on scalars: 

L Q = d^d , Li — d\di + d dQ, L 2 — d* 2 d 2 + d 2 d* 2 etc. . 
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□ 

Remarks. 

1) The case p = is special because Lq only consists of oIqcIq while L p with p > 
has two parts d p _id*_ 1 + d* p d p . The degree formulas actually count closed paths of 
length 2 starting at x. For p > 0, paths can also lower the dimension. For example, 
for p = 1, where we start with an edge, then there are two paths which start at the 
edge, chose a vertex and then get back to the edge. This explains the correction 
p+l. 

2) For p = 1, we can count the number of triangles attached to an edge x with 
deg 1 (x) = Li(x,x) + 2. 

3) A closed path interpretation of the diagonal entries will generalize the statement. 
Paths of length 2 count adjacent simplices in Q. 

We can read off the total number v p of p-simplices in G from the trace of L p . This 
generalizes the Euler handshaking result that the sum of the degrees of a finite 
simple graph is twice the number of edges: 

Corollary 14 (Handshaking). tr(L p ) = (p + 2)v p+ i. 

Proof. Summing up deg p (x) = L p (x, x)-(p+ 1) gives v p+1 = tr(L p ) - (p+ l)v P +i- 

□ 

We write just 1 for the identity matrix. The formula 

str(L + 1) = X (G) 

follows from McKean-Singer because str(l) = x(G),str(L) = 0. Combinatorically 
it is equivalent to a statement which we can verify directly and manifests in can- 
cellations of traces: 

tr(L + 1) = 2vi + v 
tr(Lx + 1) = 3v 2 + 3v 1 
tr(L 2 + 1) = 4w 3 + 4v 2 

Adding this up gives str(L+l) = tr(L + l)— tr(Li + l) + - • • = v — vi+v 2 — v 3 + - • • = 
X(G). 

An other consequence is 

tr(L )/2 - wi 
tr(Li - 2)/3 = v 2 
tr(L 2 - 3)/4 = v 3 

Remarks. 

1) Each identity str(L k + l) = x(G) produces some "curvatures" on the super graph 
Q satisfying 

X (G) = ^) 

the case k = being trivial giving k(x) = ^—i^ dlm ( x ) for which Gauss-Bonnet is 
the definition of the Euler characteristic and where k = 1 is the case just discussed. 
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While the adjacency matrix A of a graph has tr(A k ) as the number of closed paths 
in G of length k, the interpretation of tr(L k ) = tr((B — A) k ) becomes only obvious 
when looking at it in more generally when looking the full Laplacian L as we do 
here. Instead of finding an interpretation where we add loops to the vertices, it is 
more natural to look at paths in Q: 

Definition 7. A path in Q = [J Q\. as a sequence of simplices x k <G Q such that 
either x k is either strictly contained in x k+ i or that x fe+1 is strictly contained in x k 
and a path starting in Q k can additionally to Qk only visit one of the neighboring 
spaces Gk+i or Qk~i along the entire trajectory. 

The fact that this random walk on Q can not visit three different dimension-sectors 

Gk-i,Gk, Gk+i 

is a consequence of the identities d 2 — (d*) 2 = 0; a path visiting three different sec- 
tors would have cases, where d or d* appears in a pair in the expansion of (d+d*) 2k . 
Algebraically it manifests itself in the fact that L n is always is a block matrix for 
which each block L£ leaves the subspace Cl k of k- forms invariant. 

Examples. 

1) For a graph without triangles, a path starting at a vertex v is a sequence 
Vq, ei, V\, e2, i>2 • • • • Every path in Q of length 2n corresponds to a path of length n 
in G. 

2) For a triangular graph G, there are three closed paths of length 2 starting an 
edge e = (v\, i^)- The first path is e, v\, e, the second e, V2, e and the third is e, t, e 
where t is the triangle. 

3) Again for the triangle, there are 6 closed paths of length 4 starting at a vertex: 
four paths crossing two edges v,ei,v,ej,v and two paths crossing the same edge 
twice v,ei,vi,ei,v. There are 9 paths of length 4 starting at an edge. There are 
four paths of the form e, Vj,e, v k , e and four paths of the form e, v i} ej,Vi,e and one 
path e, t, e, t, e. 

Proposition 15 (Random walk in Q). The integer D^ y is the number of paths of 
length k in Q starting at a simplex x and ending at a simplex y. The trace tr(D 2k ) 
is the total number of closed paths in Q which have length 2k. 

Proof. We expand (d + d*) k . For odd k we have expressions of the form dd*d- ■ ■ dd* 
or d*d - ■ ■ d*d. For even k, we have expressions of the form dd* ■ ■ ■ d* or d*d- ■ ■ d. 
The second statement follows from summing over x. □ 

As a consequence of the McKean-Singer theorem we have the following corollary 
which is a priori not so obvious because we do not assume any symmetry for the 
graph. 

Corollary 16. Let G be an arbitrary finite simple graph. The number of closed 
paths in Q of length 2k starting at even dimensional simplices is equal to the number 
of closed paths of length 2k starting at odd dimensional simplices. 

Proof. str(D 2k ) = tr(L k |ft b ) - tr(L k |tt f ) = 0. □ 

For example, on a triangle, there are 6 closed paths of length 4 starting at a vertex 
and 9 closed paths of length 4 starting at an edge and 9 closed paths of length 4 
starting at a triangle. There are 3-6 + 9 = 27 paths starting at an even dimensional 
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simplex (vertex or triangle) and 27 = 3 • 9 paths starting at an odd dimensional 
simplex (edge). 

7. Curvature 

Finally, we want to write the curvature K(x) of a graph using the operator D. For 
a vertex x € V, denote by Vk(x) the number of K^+i graphs in the unit sphere 
S{x). The curvature at a vertex x is defined as 

fe=i 

It satisfies the Gauss-Bonnet theorem [17] 

J2k(x) = x (g), 

an identity which holds for any finite simple graph. The result becomes more 
interesting and deeper, if more gometric structure on the graph is assumed. For 
example, for geometric graphs, where the unit spheres are discrete spheres of fixed 
dimension with Euler characteristic like in the continuum, then K(x) = for odd 
dimensional graphs. This result [20] relies on discrete integral geometric methods 
and in particular on [5T] which assures that curvature is the expectation of the 
index of functions. 

Definition 8. Denote by L p (x) the operator L p restricted to the unit sphere S(x). 
It can be thought of as an analogue of a signature for differential operators. Define 
the linear operator A p — > A' p (x) = A p -i(x)/(p+l) so that A p (x) — A p ^2{x) / {p{p+ 
1)). 

Corollary 17. The curvature K{x) at a vertex x of a graph satisfies 

K(x) = str(L(x)") . 
Proof. From tr(L p (x)) = (p + 2)V p+ i we get V p -x{x) = tr(L p _ 2 (x))/p and so 

K (x) = yY^ = Y t i^A = str ( L "(x)) . 

□ 

We have now also in the discrete an operator theoretical interpretation of Gauss- 
Bonnet-Chern theorem in the same way as in the continuum, where the Atiyah- 
Singer index theorem provides this interpretation. 

If we restrict D to the even part D : fit, — > flf, the Euler characteristic is the ana- 
lytic index ker(D) — ker(D*). The topological index of D is J^xex str(D 2 (x)"). 
As in the continuum, the discrete Gauss-Bonnet theorem is an example of an index 
theorem. 
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8. DlRAC ISOSPECTRAL GRAPHS 

In this section we describe a general way to get Dirac isospectral graphs and give 
example of an isospectral pairs. In the next section, among the examples, an 
other example is given. There is an analogue quest in the continuum for isospectral 
nonisometric metrics for all differential forms. As mentioned in j3], there are various 
examples known. Milnor's examples with flat tori lifts to isospectrality on forms 
for example. 



Figure 2. The Halbeisen-Hungerbiihler pair is an example of a 
pair of isospectral graphs for the Dirac operator D. Both graphs 
have vq = 16 vertices, v\ = 21 edges and vi — 2 triangles so that 
x(G) — vq — vi + V2 = —3. The Betti numbers are bo — 1, foi = 
4, &2 = 0. Since the spectrum of Li is {3,3} for both graphs and 
Sunada methods in |llj have shown that the Laplacians Lo(Gi) on 
functions are isospectral, the McKean-Singcr theorem in the proof 
of Proposition [18] implies that also Li are isospectral so that the 
two graphs are isospectral for D. 



Proposition 18. Two connected finite simple graphs G±,G2 which have the same 
number of vertices and edges and which are triangle free and isospectral with respect 
to Lq are isospectral with respect to D. 

Proof. Since x(Gi) = vq — V\ = bo — b\ = 1 — b\ are the same, Hodge theory shows 
that L\ has b\ zero eigenvalues in both cases. Because the graphs are connected, 
I/o have bo = 1 zero eigenvalues. The McKean-Singer theorem implies that the 
nonzero eigenvalues of L match the nonzero eigenvalues of L\. □ 



This can be applied also in situations, where triangles are present and where the 
spectrum on L2 is trivially equivalent. An example has been found in jllj . That 
paper uses Sunada's techniques to construct isospectral simple graphs. We show 
that such examples can also be isospectral with respect to the Dirac operator. 
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Figure 3. An example of a cospectral pair for the Laplace op- 
erator L given by Haemers-Spence [TD]. Both graphs have v = 8 
vertices, V\ = 10 edges and v-i — 2 triangles so that x(G) = 
— v\ + v-2 = 0. The Bctti numbers are bo = 1, &i = 1,&2 = 0. 
While the spectra of Lq agree, the spectra of L2 are not the same. 
For the left graph, the eigenvalues are 4, 2, for the right graph, the 
eigenvalues are 3 and 3. We would not even have to compute the 
eigenvalues because tr(L|), which has a combinatorial interpreta- 
tion as paths of length 4 in Q does not agree. We conclude from 
McKean-Singer (without having to compute the eigenvalues) that 
also the eigenvalues of L\ do not agree. 



These examples are analogues to examples of manifolds given by Gornet [3 where 
isospectrality for functions and differential forms can be different. 




Figure 4. Two isospectral graphs for Lq found in [57] • Both are 
contractible and have Euler characteristic 1. Since they have the 
same operator L%, by McKean Singer also the spectra of L\ are 
the same. The two graphs are isospectral with respect to D. 



The rest of the article consists of examples. 



the mckean-singer formula in graph theory 
9. Examples 

9.1. Cyclic graph. Let G be the cyclic graph with 4 elements. 
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Figure 5. 



The exterior derivative d is a map from £l — R 4 — > J7 1 = i? 4 . It defines i = <^*^- 
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The bosonic eigenvalues are {4, 2, 2, }, the fermionic eigenvalues are {4, 2, 2, }. 
The example generalizes to any cyclic graph C n . The operators Lo = L\ are Jacobi 
matrices 2 — r — r* which Fourier diagonalizes to the diagonal matrix with entries 
Xk = 2 — 2cos(27rfc/n), k, 1, . . . , n. 



20 



OLIVER KNILL 



9.2. Triangle. Let G be the triangle. 



Vi=3 Vp3 Vi=l 



1 dim= 2. 




bp 1 bp bp 




Figure 6. 



The exterior derivative d : £7q — ► f^i is 
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The later has the bosonic eigenvalues {3, 3, 0, 3 } and the fermionic eigenvalues 
{3,3,3 }. We have 



tr(e~* Lo ) = 
and so str(e tL ) = 1. 



-3t I Jit 



e M + 2) ,tr(e 
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-3*. 
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9.3. Tetrahedron. For the tetrahedron K4, the exterior derivative do is a map 
from 0, = i? 4 ->■ fl 1 = R 6 . 
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do do 



and Lq = do do 
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which sums up to L± — d\d\ + dod^ = Aid. 



Finally, d 2 = [1, 1, 1, 1] and L 2 — d* 2 d 2 + d\d\ — 



For a complete 



graph K n , the eigenvalues of Lq are {0,n, . . . ,n }. The other operators Lk with 
k > 1 are diagonal with entries n. The Dirac operator D restricted to p-forms is a 
B(n + l,p+l) x B(n + l,p) matrix, where B(n, k) = n\/(k\(n— k)l) is the binomial 
coefficient. 
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9.4. Octahedron. Lets just write out the matrix L = D 2 
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Figure 8. 
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tr(-tLo) = e- 6t (3e 2 



+ e M + 2 



, e 2t 

tr(-tLi) = 3e-«(e 2t + l)" 

tr(-tLa) = e- 6t (e 2 * + l) 3 
rives the super trace 

str(-tL) = tr(- 



<t(Lo) 
ct(L 2 ) 



{0,4,4,4,6,6} 

{2,2,2,4,4,4,4,4,4,6,6,6} 

{0,2,2,2,4,4,4,6} 



•tLo) -tr(-tLi) +tr(-tL 2 ) 
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Figure 9. 

9.5. Icosahedron. We write X k if the eigenvalue A appears with multiplicity k. 
The bosonic eigenvalues of L = D 2 are 

2 , 2 5 , 3 4 , 5 4 , 6 5 , (3 - V5f, (5 - V^) 3 , (3 + V5) 3 , (5 + V5) 3 . 

The fermionic eigenvalues of L are 

2 5 , 3 4 , 5 4 , 6 5 , (3 - Vb) 3 , (5 - V5) 3 , (3 + V5) 3 , (5 + ^5) 3 . 

The eigenvalues appear only on bosonic parts matching that the Betti vector 
ft = (1,0,1) has its support on the bosonic part. The picture colors the ver- 
tices, edges and triangles according to the ground states, the eigenvectors to the 
lowest eigenvalues of Lq,Li and Li. 




Figure 10. The matrices Re(exp(iD)) = cos(£>) and 
Im(exp(z_D)) — sin(£>) which appear in the solution of the wave 
equation on the graph. Having the matrix D in the computer 
makes it easy to watch the wave evolution on a graph. 
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9.6. Fullerene. The figure shows a fullerene type graph based on an icosahedron. 
It has Vq = 42 vertices, v\ = 120 edges and Vi = 80 triangles. The Dirac operator 
D is a 242 x 242 matrix 



X= 2 dim= 2. 



b r 1 b 1= 1 





, .. ■ ■'■'..^tv- r tl .>. 
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Figure 11. 

In the next picture we see the three sorted spectra of L$, L\ and Li. 



(92,270. 180] 
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Figure 12. 
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9.7. Torus. Here is the spectrum of a two dimensional discrete torus 




16 



32 



x= o 



dim= 2. 



bri= 1 



b,= 2 



1 




and a discrete Klein bottle: 



Figure 13. 
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Figure 14. 
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9.8. Dunce hat and Petersen. The dunce hat is an example important in ho- 
motopy theory. It is a non-contractible graph which is homotopic to a one point 
graph. 




v,= 36 
2. 




Figure 15. 



Petersen graph 





Figure 16. 
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9.9. Isospectral graphs. As mentioned in the previous section, some known isospec- 
tral graphs for the graph Laplacian are also isospectral for the Dirac operator. 




Figure 17. 

In [57], all isospectral connected graphs up to order vo = 7 have been computed. 

For order vq = 7, there are already examples which are also isospectral for all p- 

forms and for the Dirac operator. Some of them lead to Dirac isospectral graphs: 

(t(D) = {-y/3 + y/2, \/3+ ^2,-^2 + V3, + y/Z, -y/S, - V?> 

\/3, y/3, -y/l -y/2,y/3- y/2, -1, 1,-^2- \/3, \/2 - y/Z, 0}. 

The Laplacian eigenvalues are a(L ) = {3 + V2, 2 + v3, 3, 3 — v2, 1, 2 — V3, 0} 

and o-(Li) = {3 + y/2, 2 + s/3, 3,3, 3 - ^2, 1,2- V^} as well as a(L 2 ) = {3}. 
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9.10. Planar regions. One of the questions which fueled research on spectral the- 
ory was whether one "can hear the shape of a drum" . The discrete analogue 
question is whether the Dirac operator can hear the shape of a discrete hexagonal 
region. 



[110,278, 168} 




Figure 18. 

A hexagonal region and the spectra of the Laplacians Lq, L±, 



<47, 116,70] 




Figure 19. 

A convex hexagonal region and the spectra of L , L\, L 2 - 

As far as we know, no isospectral completely two dimensional convex graphs are 
known. 
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9.11. Benzenoid. Isospectral benzenoid graphs were constructed in pQ. They are 
isospectral with respect to the adjacency matrix A. but not with respect to the 
Laplacian L = B — A. After triangulating the hexagons, the hrst pair of graphs have 
Euler characteristic 0, the second Euler characteristic 1. But after triangularization, 
the graphs are not isospectral with respect to A, nor with respect to L. 




Figure 20. Two examples of isospectral Benzenoid pairs which 
are isospectral for the adjacency matrix A. These graphs are one 
dimensional since they contain no triangles. McKean-Singer shows 
that the nonzero spectrum of Lq is the same than the nonzero 
spectrum of L\. 




78 Vf 102 




FIGURE 21. The Dirac spectrum of the first Benzenoid. 
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9.12. Gordon- Webb- Wolpert. Isospectral domains in the plane were constructed 
with Sundada methods by Gordon, Webb and Wolpert [5] . A simplified version con- 
sists of two squares and three triangles. 

Figure 22. 




Figure 23. 

To try out graph versions of this, we computed the spectra of a couple of graphs of 
this type. They are all homotopic to the figure 8 graph but are all not isospectral. 
This shows that discretizing Sunada can not be done naivly. The next section men- 
tions examples given in [8] which were obtained by adapting Sunada type methods 
to graphs. 
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9.13. Halbeisen-Hungerbiihler. Here are the picture of the Dirac matrices and 
their spectrum of the isospectral graphs constructed in [8]. 




Figure 24. 
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9.14. Polyhedra. For the 5 platonic solids, we have the following Dirac spectra 
and complexities JIa^o ^> wnere we write A™ for A with multiplicity n. 



Solid 


Non-negative Dirac Spectrum 


Complexity 


Tetrahedron 




_ 2 14 


Octahedron 




2 18 • 3 3 


Cube 


U/e^v^o 6 } 

{v^v^v^Vs-v^V} 


-2 10 • 3 


Dodecahedron 


-2 11 • 3 4 • 5 4 


Icosahedron 


{v / 2 5 ,v / 3 4 ,v / 5 4 ,v / 6 5 ,\/5±v / 5",v / 3±v / 5 3 ,0 2 } 


2 22 • 3 9 • 5 7 



The following table confirms that for two-dimensional geometric graph for which 
the unit sphere is a circular graph, the Dirac complexity is positive. We actually 
came up with the statement after watching the following tables: 



" dim(G) 


complexity 


X(G) 


2 


-2.6121 x 10 iu 


-4 


1 


-2.6147 x 10 42 


-60 


1 


-5.9608 x 10 17 


-24 


2 


-2.1859 x 10 25 


-10 


3/2 


-4.2166 x 10 40 


-40 


3/2 


-4.7406 x 10 16 


-16 


2 


-2.1456 x 10 28 


-4 


2 


-8.2030 x 10 69 


-10 


5/3 


-5.1018 x 10 12 


-4 


5/3 


-1.0423 x 10 30 


-10 


1 


-2.2518 x 10 22 


-30 


1 


-2.4277 x 10 9 


-12 


5 

L 3 


-5.8320 x 10 6 


-2 



,C- 



dim(G) 



complexity x(G) 



1 


-4.6449 x 10 b 


-10 


2 


3.5323 


X 


10 84 


2 


2 


1.9246 


X 


1Q 35 


2 


1 


-6.6871 


X 


10 15 


-28 


1 


-1.2496 


X 


10 31 


-58 


1 


-7.8275 


X 


10 12 


-22 


1 


-3.4164 


X 


10 15 


-22 


1 


-4.0315 


X 


10 37 


-58 


3 


3.7873 


X 


10 26 


2 


3 


1.4404 


X 


10 64 


2 


2 


2.7042 


X 


10 44 


2 


2 


3.4381 


X 


10 18 


2 


3 


-1.6987 


X 


10 14 


1 




Figure 25. The DisdyakisTriacontahedron is the Catalan solid 
with maximal complexity. 
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Appendix: The original proof 

While the super symmetry argument is elegant, the original proof is illustrative 
and shows better the link to Hodge theory. We follow [21] but change notation and 
expand some of the steps slightly. Denote by E p the eigenspace of the eigenvalue 
A of L on QP{G). 

Proposition 19 (McKean-Singer). ^(-l^dim^) = f or allp> 0. 

Proof. If A is a positive eigenvalue, we just write E p instead of E p x . 

(i) EP = dEP- 1 ® d*EP +1 

Proof. That the right hand side is included in the left hand side follows from 
[d, L] = [d* , L] = 0. For example: if / = dg with Lg — Xg, then Lf = Ldg = dLg — 
dXg = Xdg = Xf and / € EP. 

Given / € dEp- 1 ^ € d* EP+ 1 , then (f,g) = (dh,d*k) = (d 2 h,k) = so that the 
two linear spaces on the right hand side are perpendicular. Assume we are given / G 
E p which is perpendicular to both subspaces. Then = (/, c£E p_1 ) = (d*f,E p ^ 1 ) 
and = (/, d*EP +1 ) = (df, EP+ 1 ) so that d* f = df = and Lf = which together 
with / e E p implies / = 0. 

(ii) Summing the identity dim(E p ) = dim(dE p_1 ) + dim(d*E p+1 ) obtained in (i) 
gi ves E P (- 1 ) Pdim ( EP ) = 

= dim(E p ) - (dimtdEP- 1 ) + dim(d*E p+1 )) 

p even podd 

= dim(E p ) - dim(dE p ) - dim(d*E p ) 

p even 

= dim(dE p - 1 ) + dim(d*E p+1 )-dim(dE p )-dim(d*E p ) 

p even 

= J2 dim(dE p - 1 ) + dim(d*E p+1 )-dim(dd*E p+1 )-dim(d*dE p - 1 ) 

p even 

= [dim(dEP- x ) - dim(dd*E p+1 )] + [dim(d*E p+1 ) - dun(d*dE p - 1 )] > . 

p even 

(iii) E P (-l) P dim(E p ) < 

From (ii) , we recycle J2 P ( ~ 1 ) p dim(E p ) = ^ p cvo „ dim(E p ) - dim(dE p ) - dim (d* E p ) . 
The claim follows now because each term is < 0: 

dim(E p ) - dim(dE p ) - dim(d*E p ) 

= dim(LE p ) - dim(dE p ) - dim(d*E p ) 

< dim(dd*E p ) + dim(d*dE p ) - dim(dE p ) - dim(d*E p ) 

= [dim(dd*E p ) - dim(dE p )] + [dim(d*dE p ) - dim(d*E p )] < . 

The two parts (ii) and (iii) imply the proposition. □ 

Corollary 20. str(e~ tL ) = x(G). 



Proof. Proposition (19) implies str(L k ) = for k > 0. A Taylor expansion of the 



super trace of the heat kernel using the proposition gives 

str(e- tL ) = str(l - t^ + t 2 ^ - • • • ) = str(l) = X (G) 



□ 
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Appendix: Hodge theory for graphs 

For t — we have by definition the Euler characteristic and for t — > oo, the 
supertrace of the identity on harmonic forms. The proof of the Hodge theorem 
stating that H P (G) is isomorphic to the space Eq(G) of harmonic p- forms is clearly 
visible in McKean-Singers proof: 

Lemma 21. a) d*L = Ld* — d*dd* and dL = Ld = dd*d. 

b) Lf = is equivalent to df — and d* f = 0. 

c) f = dg and h = d*k are orthogonal. im(d) U im(cT) span im(L). 

Proof, a) is clear from d 2 = (d*) 2 = 0. 

b) if Lf = then = (/,£/) = (d*f,d f ) + (df,df) shows that both d* f = and 
df — 0. The other direction is clear. 

c) (/, h) = (dg, d*k) = (ddg, k) = (g, d*d*k) =0. □ 
Any p-form can be written as a sum of an exact, a coexact and a harmonic p-form. 
Lemma 22 (Hodge decomposition). There is an orthogonal decomposition 

= im(L) + ker(L) = im(rf) + im(cP) + ker(L) . 

Any g can be written as g — df + d*h + k where k is harmonic. 

Proof. The operator L : — > A is symmetric so that image and kernel are per- 
pendicular. We have seen before that the image of L splits into two orthogonal 
components im(c?) and im(d*). □ 

Theorem 23 (Hodge- Weyl). The dimension of the vector space of harmonic k- 
forms on a simple graph is b^. Every cohomology class has a unique harmonic 
representative. 

Proof. If Lf = 0, then df = and so Lf = d*df = 0. Given a closed n-form 
g then dg = and the Hodge decomposition shows g = df + k so that g differs 
from a harmonic form only by df and so that this harmonic form is in the same 
cohomology class than /. We can assign so to a cohomology class a harmonic form 
and this map is an isomorphism. □ 

Appendix: Cech cohomology for graphs 

We have mentioned that the Dirac operator has appeared when computing Cech 
cohomology of a manifold and investigating the relation between the spectrum of 
the graph and the spectrum of the manifold [53]. We want to show here that 
Cech cohomology is equivalent to graph cohomology. This is practical: while graph 
cohomology is easy to implement in a computer, it can be tedious to compute 
the operator D and so cohomology groups. If the zero eigenvalues and so the 
cohomology is the only interest in D, then one can proceed in two different but 
related ways: One possibilities to reduce the complexity of computing the kernel of 
D is to deform the graph using homotopy steps to get a simpler graph. Another 
possibility is to look at the Dirac operator of the nerve. The following definition is 
equivalent to a definition of Ivashchenko: 

Definition 9. A graph is called contractible in itself, if there is a proper sub- 
graph H which is contractible and which is the unit sphere of a vertex. The graph 
G without the vertex x and all connections to x is called a contraction of G. For a 
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contractible graph, there is a sequence of such contraction steps which reduces G 
to a one point graph. 

This allows now to consider Cech cohomology for graphs: 

Definition 10. A Cech cover of G is a finite set of subgraphs Uj which are con- 
tractible and for which any finite intersection of such subgraphs is either empty or 
contractible. A Chech cover defines a new graph N called the nerve of the cover. 
The vertices of N are the elements of the cover and there is an edge (a, b) if both the 
subgraphs a, b have a nonempty intersection. Cech cohomology is then defined 
as the graph cohomology of the nerve of G. 

Remark. The number of vertices of the nerve of G is an upper bound for the 
geometric category gcat(G) of G and so to strong category Cat(G) [T5] . 



Proposition 24. The nerve of a Cech cover of G is homotopic to G. 

Proof. In a first step we expand each of the Uj so that in each Uj , the intersections 
Uk H Uj are in a neighborhood B(xj) of a point Xj € Uj. Now homotopically shrink 
each Uj to B{xj). □ 

Corollary 25. Cech cohomology on a finite simple graph is equivalent to graph 
cohomology on that graph. 

Proof. Graph cohomology obviously is identical to Cech cohomology if we chose 
the cover Uj — B(xj) for the vertices Xj which has the property that the nerve of G 
is equal to G. Because the nerve N of G is homotopic to G and graph cohomology 
of homotopic graphs is the same (as already proven by Ivashchenko |12j). we are 
done. □ 
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